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Aero-engine structures have very low inherent damping and so artificial 
damping is often introduced by pumping oil into annular gaps between the 
casings and the outer races of some or all of the rolling-element bearings 
supporting the rotors. 

The thin oil films so formed are called squeeze film dampers and they can 
be beneficial in reducing rotor vibration due to unbalance and keeping to 
reasonable limits the forces transmitted to the engine casing. However, 
squeeze-film dampers are notoriously non-linear and as a result can introduce 
into the assembly such phenomena as subharmonic oscillations, jumps and combi- 
nation frequencies. 

The purpose of the research described in this paper is to investigate such 
phenomena both theoretically and experimentally on a test facility reproducing 
the essential features of a medium-size aero engine. 

The forerunner of this work has been published in refs. [1],[2]. It was 
concerned with the examination of a squeeze-film damper in series with housing 
flexibility when supporting a rotor. The structure represented to a limited 
extent the essentials of the projected Rolls Royce RB401 engine. That research 
demonstrated the ability to calculate the oil-film forces arising from the 
squeeze film from known motions of the bearing components and showed that the 
dynamics of a shaft fitted with a squeeze film bearing can be predicted rea- 
sonably accurately. 

An aero-engine will normally have at least two shafts and so in addition 
to the excitation forces which are synchronous with the rotation of one shaft, 
there will also be forces at other frequencies from other shafts operating on 
the squeeze-film damper. The present paper is concerned with theoretical and 
experimental work to consider severe loading of squeeze-film dampers and to 
include these additional effects. The theoretical methods are similar to those 
discussed in references [1] and [2]. 


* Now with General Electric Company, Cincinnati, U.S.A. 
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NOTATION 


a 
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c 

d 

D 

F 

k 

4 

L 


m e Cw 2 

m 


m el ,2 
m h 


p cav 

p s 

P cl,2 

R 

r l,2 

T 


x 

U 1 , 2 
W 


Mass ratio m e 2/ m h 
Linear damping coefficient 

Radial clearance of the squeeze-film bearing 
Seal end-gap 

Diameter of squeeze-film bearing 

Amplitude of force transmitted to engine frame 

Stiffness of the flexible bars 6 in Fig. 2 

Squeeze-film bearing land width 

Radial depth of sealing ring 

Width of groove 

Normalising parameter for forces 
Rotor mass 

Effective mass of rotor 1,2. 

Effective mass of rotor 1 plus bell housing 4, 
squeeze-film bearing housing 5 and flexible bars 6 in 
Fig. 2. 

Cavitation pressure 
Oil supply pressure 
Unbalance force on rotor 1, 2 
Radius of the squeeze-film bearing 

Radius of unbalance of mass m^ 0 from axis of rotation. 

Transmissibility = force reaching ground unbalance 
force 

Amplitude of rotor vibration 

Effective radius of unbalance of rotor (1,2) mass 
from axis of rotation 

g/cw n 2 
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wj / <^2 


X 

♦ 

c 

M 

Hu n (R/c) 2 
W 1 , 2 



m e2 w n c 

Boundary-condition factor 
Attitude angle 
Eccentricity ratio 
Dynamic viscosity 

Normalising parameter for pressures 
Angular velocity of rotor 1,2. 
Undamped natural frequency 
Differentiation with respect to wt 


THE TEST FACILITY 

Fig. 1 shows diagrammatical ly a section through the front end of a modern 
aero-engine, having two nested rotors. The structure is subjected to a variety 
of forces, for example forces arising from unbalance in the H.P. rotor are 
transmitted into bearing housings supporting the L.P. rotor, which often incor- 
porate squeeze film dampers. A test facility was designed to reproduce such 
essential features of the aero-engine assembly, in particular to examine the 
vibration of a rotor employing a squeeze-film damper, whose housing can be 
under the influence of an excitation force which is different in frequency from 
that of its shaft rotation. 

A drawing of the test facility is shown in Fig. 2. The assembly consists 
of two sub-assemblies comprising rotors 1 and 2, representing the aero engine 
HP and LP shafts respectively. These are capable of being driven at different 
speeds by variable-speed motors through pin-and-cord couplings. The unbalance 
mass 3 produces a dynamic force which is transmitted through the bell-housing 4 
to the housing of the squeeze-film damper 5, which is itself flexibly supported 
at 6 and contains the roller bearing 7. The outer race of this bearing is 
surrounded by the squeeze-film damper 3, which is provided with either RIO oil 
or Shell calibration fluid. These have viscosities at room temperature of 15cp 
and 5cp respectively. The unbalance mass 9 provides a further force at the 
inner race of the roller bearing, whose movement is permitted by the self- 
-aligning bearing 10 and whose location represents a node of rotor 2. 

The damper has two bearing lands which are separated by a central circum- 
ferential oil-supply groove. End-plates are attached at both ends of the 
squeeze film to afford some sealing. The axial end clearance between the 
plates and the dampers and thus the extent of the sealing can be varied by the 
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insertion of spacing shins. The rotation of the damper journal is prevented by 
dogs located in one of the end plates, providing sufficient clearance to allow 
the journal to move in any part of the squeeze-film clearance space. The 
squeeze-film damper diameter is set by the outer race of the roller bearing 7 
(140 mm). The damper land length and its radial clearance are typical of aero 
engine application (Fig. 3). 


The rig parameters 

Preliminary impulse tests were carried out to define the total effective 
stiffness of the flexible bars 6 and to find the effective masses of rotors 1 
and 2 at the squeeze-film location. From tests on the flexibly-mounted housing 
alone, its natural frequency was found to be 192 Hz. Using an effective 
housing mass of 3.9 kg the total stiffness of the flexible bars was found to be 
about 5.68 x lCr N/m. With rotor 2 alone in position, the natural frequency 
was 75 Hz in the horizontal and 71 Hz in the vertical direction. The differ- 
ence was almost certainly due to the movement of bearing 7 within its squeeze 
film annulus. This was empty of oil at the time of testing and promoted bounc- 
ing of bearing 7, By taking 73 Hz as an averaged frequency, the effective mass 
of rotor 2 was found to be 22.16 kg. A similar test was carried out with rotor 
1 in position and with rotor 2 removed, and finally with both rotors in posi- 
tion. Table 1 shows the findings of all the tests in a concise form. 


Table 1 


1 

1 

! Arrangement 

i 

i 

i 

i 

i 

i 

Horizontal 

(Hz) 

Vertical 

(Hz) 

Effective 

mass 

(kg) 

Accepted 
Natural 
Freq. ; 
(Hz) ; 

i 

i 

! Squeeze film 
! (SF) housing 
! only 

i 

192 

192 

3.9 

192 : 

i 

! Rotor 2 and 
! SF assembly 

t 

i 

i 

75 

71 

22 . 2 

(rotor 2 
only) 

73 ; 

i 

! Rotor 1 and 
1 SF assembly 

t 

i 

1 

94 

103 

10.0 

(rotor 1 
only) 

98 , 

i 

! Rotors 1 and 2 

l 

i 

61.2 

63.6 

32.4 

62 ; 


The undamped natural frequency of the complete assembly was taken as 62 Hz 
(3720 c/min) and it can be seen that the effective mass in the last case is 
about equal to the sum of those in the previous two cases. The static load 
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provided by rotor 2 on the squeeze-film was half its weight, corresponding to 
50/2 or 25 kg. The dynamic loads arose from the unbalance masses 3 and 9, and 
can be expressed as P cl = mjrjuj" and P c 2 = m^r-,!^". where mj and m 2 are re- 
spectively the unbalance masses on rotors 1 and 2, rj, are the radii of the 
mass locations from the centre of rotation of each shaft, and and are the 
speeds of the shafts. However, in the analysis it is more convenient to con- 
sider the unbalance mass-radius product to have a magnitude equal to the total 
effective mass of each rotor at an effective unbalance radius u. 


SINGLE-SHAFT VERSION OF THE TEST FACILITY 

The test facility was first separated into its two rotor components by 
removing the bell housing 4. Tests were carried out on the LP rotor (2) alone 
mounted in the squeeze-film damper to investigate its performance characteris- 
tics . 


A diagrammatic representation of the test facility in this configuration 
is given in Fig. 4. The degree of damping offered by the squeeze-film is quite 
critical. If it is too great the damper operates rather like a rigid link and 
conveys little benefit to the system. On the other hand if the damping is too 
small, excessive movements take place in the squeeze film annulus. 

Assuming for a moment that the damping is linear, its effect in series 
combination with the rotor mass and housing flexibility is shown in Figs. 5a, 
b, c. From these figures it can be seen that there is likely to be an optimum 
value of damping b to ensure reasonable levels of rotor vibration relative to 

the engine frame (Fig. 5a), relative to the bearing housing (Fig. 5b) and of 

transmissibility (Fig. 5c). Complications are introduced by the fact that due 
to non-linearity of the squeeze-film its damping increases greatly with eccen- 
tricity ratio e. Fig. 5b indicates that high values of e can arise over a 

considerable frequency range, particularly at low values of b. Under these 
conditions the apparent benefits of low damping indicated in Figs. 5a and 5c 
may not be realised. 

The most critical operating condition occurs at a speed corresponding to 
the natural frequency of the assembly, “2 is 1 . From Figs. 5a,b,c, three dis- 
tinctive features can be observed. "n 

a) The non-dimensional amplitude of vibration of the rotor x 0 /u 
increases with non-dimensional damping b/Vkiri. 

b) The non-dimensional eccentricity of the rotor in the squeeze- 
film damper, e/(u/c), is independent of non-dimensional damping 
b//km, where e = (x 0 -Xj)/c. 

c) The transmissibility T increases with non-dimensional damping 
b/./km. 

Thus the assembly is quite unlike one having a damper in parallel with its 
flexibility . 
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Preliminary tests were first run on the reduced test facility to ascertain 
whether these trends were exhibited in the actual non-linear system. Various 
means were adopted by which the effective damping could be altered. For 
example, damping could be increased by 

i) increasing the oil viscosity 

ii> tightening the sealing, or 

iii) increasing the supply pressure. 

Figs. 6a, b, c show the results obtained for various sealing conditions, 
oil viscosities and supply pressures. From these, conclusions can be drawn 
which are consistent with the three distinctive features using the linear 
damping model. This model thus provides a reasonable qualitative guide to 
damper performance in the most critical region of rotor operation. 

The next purpose of the reduced test facility was to look for evidence of 
jump phenomena and sub-harmonic vibrations in the frequency responses as 
these had been observed in certain engine tests where the oil viscosity was 
low. 


Efforts were therefore made to produce conditions on the test facility, 
which were likely to be conducive to such jumps. This indeed proved to be 
possible when the support flexibility was locked out and when the supply 
pressure to the squeeze-film was at a fairly low level. Fig. 7a shows a 
sequence of orbits of rotor vibration relative to the squeeze-film damper 
housing using a supply pressure of 5 psi. At about 3450 rev/min a pronounced 
jump up was observed. This was almost certainly due to a change in the 
effective parameters governing system performance. This change was probably 
brought about by sudden venting of the squeeze film clearance space from atmos- 
phere, resulting in a change of level at which the negative hydrodynamic 
pressure in the squeeze film was curtailed. Fig. 7a shows a set of pressure 
recordings which seem to bear this out. This would result in much smaller 
squeeze-film forces and hence a much larger orbit of vibration of the rotor 
centre, and can be predicted by numerical computation (Fig. 7b). For the 

latter purpose, the hydrodynamic squeeze-film forces were numerically derived 
by integrating the oil-film pressure around the circumference of the damper. 
These forces depend on the oil-film cavitation model used and on the end-plate 
clearance expressed as a X factor [2]. The X-factor describes the boundary 
conditions imposed at the ends of the squeeze-film damper by the tightness of 
the end-seals and in the present research has been extended to cover various 
1 /q ratios (Fig. 7c). In the present case where the J l/ D ratio is (1045 a 

X-factor of 0.09 was used, being appropriate to an end plate gap/clearance 
ratio of 0.23 and an end-plate clearance of 0.053 mm. The non-linear equa- 
tions of motion were then solved by integration at each time step of a marching 
scheme. 


With regard to subharmonic vibrations. Fig. Sa shows a succession of 
orbits of the LP shaft (2) relative to the squeeze-film housing for a u/c value 
of 0.33 and using Shell calibration fluid as the lubricant. The speed range 
covered was 3330 to 5075 rev/min. with the undamped natural frequency at 4300 
rev/min (73 Hz). The centre of the LP rotor executed a butterfly type of 
motion, such that it completed each "wing" for alternate revolutions of the 


66 



rotor. This half-engine order motion persisted over a wide range of opera- 
tional speed. Eventually the "wings" coalesced after the critical speed was 
passed to produce an essentially single orbit. The predicted orbits for the 
same speeds are shown in Fig. 8b, in which the same double looping is clearly 
evident. Coalescing is predicted somewhat later in the speed range than was 
measured. 


Z -MODULATION 

The z-mod technique presents a great deal of useful information on the 
harmonic content of a signal in a very compact way and by this means it is 
possible to appreciate an overall performance signature of an engine over a 
wide speed range. A z-mod is similar to a "waterfall" diagram used in rotating 
machinery diagnostics except that engine speeds are arranged along the x 
axis, instead of along the y axis, frequency content is displayed along the 
y axis instead of along the x axis and the strength of a frequency component 
is represented as a colour (or as a degree of blackness) instead of a laid-back 
peak. 


By way of elucidating this technique, consider Fig. 3a which shows an 
experimental z-mod of horizontal components of LP rotor vibration relative to 
the squeeze-film housing using a specially-adapted commercial spectral 
analyser. It will be observed that half-integer order frequency components 
manifest themselves strongly around the natural frequency of 73 Hz. 

Software was developed to obtain theoretical z-mods with the wide range 
of colours available on many microcomputers. This made possible direct com- 
parisons of predictions with the experimental z-mods. 

The z-mod presentation of numerically-predicted horizontal LP shaft vibra- 
tion relative to the squeeze-film damper housing is shown in Fig. 9b, and 
compares well with the experimental z-mod of Fig. 9a. Half engine orders are 
indeed predicted around 73Hz. This is due to squeeze-film non-linearities 
manifesting themselves when the overall vibration becomes appreciable, as it 
does around the natural frequency of 73Hz. 

Half-integer orders are a common feature in loose assemblies [3], It is 
not unreasonable to suppose that in the present case the squeeze-film damper 
clearance was acting as a loose assembly on account of the fairly high unbal- 
ance load together with the fairly short squeeze-film damper lands and the low 
viscosity of the oil supplied. When heavier oil was introduced all half-integer 
order harmonics disappeared. 

A non-linear system can be conceived in which a subharmonic may be 
excited, whose frequency is half that of the excitation frequency. To show that 
this is possible, consider a system with a square-law restoring force subject 
to a steady force (such as gravity) and to harmonic excitation (such as that 
due to unbalance), both of which are provided by the test facility. Let the 
system satisfy the differential equation, 

.. o 

mx + yx + ftx = P Q + P| cos wt. say 
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that is, the effective stiffness of the squeeze film is y-r ux. 
Let us assume that a possible solution is 
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A cos - t. 
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This equation is satisfied for 
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The first two conditions imply that P Q = Pj . which happens to be conducive 
to rattle since then the gravity force P Q is just neutralised by the dynamic 
lift-off force Pj. This example shows that a second harmonic, resulting from 
the non-linear term, can prevail, being compensated by a certain intensity of 
excitation at some specific frequency. This frequency is illustrated most 
strikingly when a squeeze-film damper is accompanied by a parallel spring of 
stiffness y. In that case u corresponds to twice the first undamped natural 
frequency and a subharmonic resonance of order h is exhibited [4j. 


Since both harmonics and subharmonics can occur simultaneously, oscilla- 
tions having any half integer ratios to the excitation frequency are possible, 
given the correct form of non linearity. These oscillations are in fact combi- 
nation frequencies and will be discussed later. 


THEORETICAL ANALYSIS FOR COMPLETE CONFIGURATION OF THE TEST FACILITY 

The dynamical ly-equivalent structure of the test facility with both rotors 
connected is shown in Fig. 10. It may be shown [2] that a set of non-dimen- 
sional parameters which describe this non-linear system can be expressed as W, 

P n , u 1 - 2 / c ’ a ’ 01 » u 2 ^ u n ‘ p cav /yJw n (R/c) ‘'’ and P s /Hw n (R/c> 2 . 

The non-linear equations of motion were solved by again integrating at 
each time step of a marching scheme. To obtain a feel for the kind of results 
obtained. Fig. 11a shows a typical vibration orbit for the housing of the 
squeeze-film bearing relative to ground while Fig. lib shows the vibration 
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orbit of rotor 2 relative to ground. In Fig. 11c the orbit of rotor vibration 
relative to the squeeze-film housing is shown. For these cases the following 
values were applied: 

No-end sealing W = 0.2S a = 2.0 

p n = 0.016 a = 1.5 

P'jjjlR/c)* = 154.5 psi w 2^ w n = 

Uj/c = 0.23 p cav = -20 psi 

u 2 /c = 0.23 p s = 20 psi 

By putting ot = -2.0, counter-rotation of rotor 1 relative to rotor 2 is 
represented. The resulting orbits for this condition are shown in Fig. 12 in 
the same sequence as for Fig. 11. The same cases were then treated with end- 
plates included (Fig. 13). These last orbits indicate a much more stable system 
and steady states are achieved much more quickly, due to the higher damping 
ensured by the end plates. 


COMPARISON OF NUMERICAL ORBITS AND EXPERIMENTAL OBSERVATIONS FOR 
COMPLETE CONFIGURATION OF THE TEST FACILITY 

Fig. 14 shows a specimen set of comparisons of rotor 2 vibration orbits 
relative to the ground, when rotating at 3000 rev/min with rotor 1 stationary. 
For these comparisons the full non-linear equations were solved again using 
the X-f actor approach developed in ref. [2]. Fig. 14 is for a case of fairly 
tight sealing (X = 0.0356) and the agreement between experimental and pre- 
dicted orbits of vibration is seen to be reasonable, the scale factors used for 
each set of orbits being the same. This agreement is representative of com- 
parisons over a wide range of speeds, oil viscosities, supply pressures and 
unbalances and a cavitation pressure of -20 psi has been used throughout. 

Consider now cases in which rotors 1 and 2 rotate in first co-rotation and 
then counter-rotation relative to each other. Detailed comparisons between 
predicted and experimental orbits are not appropriate here due to the indeterm- 
inacy of the phasing between shafts 1 and 2 arising from slight drifts in 
speed. However it is possible to observe certain distinctive features from 
comparisons of experimental orbits of rotor vibration relative to ground for 
the co-rotating and counter-rotating cases. The same distinctive features can 
then be sought in the comparison of the numerically-predicted orbits. 

Some experimental orbits are shown in Figs. 15 and 16 (upper) for uj/u 2 = 
±1.3. The most distinctive differences are observed to be the rounded nature 
of the orbits for the co-rotating cases as opposed to the spiky nature for the 
counter-rotating. The same features can be observed in the numerically-pre- 
dicted orbits. Fig. 15 and 16 (lower). In addition, there is broad agreement on 
orbit size. 
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It was noticeable during running that the test facility was much more 
noisy for counter rotation, suggesting that transmitted forces were higher, 
perhaps due to the more sudden reversals in the orbits. 


FREQUENCY CONTENT OF SYSTEM VIBRATION WHEN SHAFTS AP.E IN CO-ROTATION 
AND IN COUNTER-ROTATION : COMBINATION FREQUENCIES 

In a linear system resonance occurs when the system is excited by a har- 
monic force of frequency near to the undamped natural frequency of the system. 
In a non-linear system other kinds of resonance are possible. 

Consider a non-linear system subject to an excitation consisting of two 
harmonic components as we have in this test facility. Thus, 

mx i- f(x) = Pj cos (wjt * 0j) + Pt cos (w 0 t + 0 9 ) 

in which we can expand f(x) in the form of a Taylor series, 

o 3 

f(x) = ajX + a 2 x" + a-jX + ••• 

Now assume that x is roughly harmonic, that is 
x = Xj cos wjt + X 2 cos w 2 t 

Thus f(x) will consist of terms in 
COS Uj -jt 


9 

COS” Wj 2 t = 

(1 

+ cos 2 J, ot) 
1 » — 


COS w ^ , o t — 

[3 

cos o ^ cos 3 ot ] 


COS Wjt COS w^t 

= 

^ COS (uj + ) t + ^ COS ( '« 

— Ulr) ) t 

COS 2u^t COS w 0 t 

= 

cos + uo't * ** cos 

( - wo)t etc. 


All possible linear combinations of the two frequencies uj and t5lus 
appear and in general it can be said that the solution will contain the fre- 
quencies 

nwj 

and nwj ± mu 2 

where m and n are whole numbers, the latter being known as combination 
oscillations. They also manifested themselves in half-integer orders during 
tests on the single-shaft version of the test facility (Fig. 9). 

Fig. 17 shows experimental z-mods for cases of co-rotation and counter- 
rotation in which wj/uo = ±1.3. These relate to vertical vibration of the 
squeeze-film housing relative to ground. Many LP and HP engine orders are 
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present in both, as are probable combination frequencies. By taking a vertical 
section of the z-mod at L.P. rotor speeds of say 40, 60 and 75 rev/s, frequency 
spectra can be constructed. Figs. ISa, b show such spectra for co-rotation 
and counter-rotation respectively. Some suggestions are made on Fig. IS as to 
what these frequencies might be in terms of HP and LP rotor speeds. 


CONCLUSIONS 

The purpose of the research described in this paper has been to reproduce 
on a test facility phenomena exhibited in aero-engine tests. These included 
subharmonic oscillations, jumps in frequency response and combination frequen- 
cies. All of these effects were reproduced and could be predicted by numerical 
computation or by appreciation of non-linear mechanics. Furthermore, at speeds 
in the vicinity of the undamped natural frequency, a linear treatment gave a 
sound understanding of trends in vibration performance. 

A method was presented by which the above features could be appreciated in 
a comprehensive way. This is known as the z-modulation technique and utilises 
up-to-date computer graphics facilities. 
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Fig. 2 : Experimental test facility. 


1 Rotor 1 
? Rotor? 

3 Unbalance moss 
U Bell housing 

5 Squeeze - film hooting 

6 Flexible bars 

7 Roller bearing 

0 Squeeze- Mm domper 
9 Unbalance mass 
X) Self -aligning bearing 



c = 0.23 mm 
d *= 0.076 ram 
D = 140 ram 
g = 0.63 mm 
i * 6.26 mm 
Lg * 7.00 mm 

L « 7.76 mm Fig. 4 - Mathematical model of the rig. 

Fig. 3 : Damper geometry. 
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a) Experimentally-observed orbits 



3080 rev/min 4059 rev/min 

b) Numerically-predicted orbits 


5075 rev/min 


Fig. 8 : Subharmonic vibrations. 
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Fig. 10 - Line diagram of test rig. 
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Fig}3: Vibration orbits for the cases of Figs. 5, 6 

but with end plates 
X* O-i. 
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Rotor vibration relative to ground 
Figl-4-.Rlo oil shaft 2 rotating at 3000 rev/min 
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Fig. 18 : Frequency spectra. 






